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ABSTRACT 

It  la  ^30\m  that  the  probability  of  receiving  light  from  an  object 
viewed  through  a  turbulent  atmosphere  follows  a  normal  Gaussian  distribution. 
Furthermore  the  root  mean  square  angular  deflection  of  the  points  of  any 
object  will  be  proportional  to  the  square  root  of  the  object-to-observer 
disteince. 

From  relations  of  the  type  described  in  the  examples,  it  is  possible  to 
predict  the  apparent  contrast  throughout  a  given  scene,  provided  the  inherent 
contrast  distribution,  the  Optical  Air  State,  and  the  range  of  the  target  is 
loiown.  The  Optical  Air  State  for  a  given  condition  of  atmosphere  can  be 
measured  using  a  telephotometer  and  a  series  of  long  thin  black  bars  of  varying 
widths. 

*This  ^per  was  prepared  under  a  contract  between  the  Universitj'"  of 
California  and  the  Bureau  of  Ships,  U.  S.  Navy.  Contribution  from 
The  Scripps  Institution  of  Oceanography,  University  of  California, 

New  Series  No.  _ . 

Present  address:  Rand  Corp.,  Santa  Monica,  California 

■***  Present  address:  Hughes  Aircraft  Co,,  Culver  City,  California 


INTRODUCTION 


One  of  several  factors  determining  the  form  of  the  image  of  a  distant 
object  is  the  time  varying  distortion  caused  by  inhomogenities  in  the 
refractive  index  of  the  atmosphere.  As  is  well  known,  the  shimmer  of  distant 
objects  is  caused  by  spatial  and  temporal  variations  in  tho  index  of  refrac¬ 
tion  of  rising  columns  of  hot  air,  mth  the  resulting  disturbance  of  the 
propagated  electromagnetic  radiation.  This  disturbance  causes  time  varying 
distortion  of  the  image  and  loss  of  fine  detail* 

This  phenomenon  is  of  particular  interest  to  the  geodesist,  because  it 

limits  the  precision  xvith  which  a  telescope  can  be  pointed  at  a  distant 
1 

object  j  and  to  the  astronomer,  because  it  blurs  the  images  of  the  stars  and 
limits  his  resolving  pov/er. 

Experiments  on  time  varying  characteristics  of  optical  transmission 
through  an  inhomogeneous  atmosphere  have  been  made  by  Riggs,  Mueller,  Graham, 
and  Mote  ,  who  found  that  for  a  path  length  of  several  hundred  meters  the 
average  apparent  instantaneous  displacement  of  a  point  on  a  target  was  greater 
than  3  seconds  of  arc  with  a  maximum  displacement  of  9  seconds.  The  instan¬ 
taneous  displacements  of  points  separated  by  more  than  5  minutes  of  arc  were 
uncorrelated. 

Some  of  the  time  varying  characteristics  of  the  shimmer  effect  are  slow 
enough  to  be  discernable  to  the  eye,  but  much  of  the  effect  is  so  rapid  that 
only  time  averaged  values  of  the  radiance  of  the  image  are  important  for 
visual  observation. 

The  analytic  treatment  of  the  propagation  of  electromagnetic  waves 

1  Washer,  F.  E,,  Williams,  H.  B.,  "Precision  of  Telescope  Pointing  for  Outdoor 
Targets,"  J,  Opt.  Soc,  Am,  36,  400  (1946), 

^  Riggs,  Mueller,  Graham,  and  Mote,  "Photographic  Measurements  of  Atmospheric 
Boil,"  J.  Opt,  Soc.  Am.  37,  415  U947). 
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through  inhomogeneous  media  has  been  considered  by  a  large  nuniber  of  investi¬ 
gators.  The  rigorous  solution  of  the  problem  by  perturbation  methods  in 
terms  of  the  statistical  parameters  of  the  inhomogenities  yields  results  in 
very  untractable  forms,  especially  for  one  interested  in  the  transmission 
of  optical  images.  Results  derived  from  approximate  models  have  been  useful 
for  specific  problems.  The  use  of  the  approximation  of  geometric  optics 
yields  results  in  a  form  that  are  most  easily  used  in  problems  of  optical  image 
transmissionj  particularly  for  the  consideration  of  the  time-averaged  image* 

The  validity  of  the  geometrical  optics  approximation  of  a  particular 
case  can  be  ascertained  by  the  value  of  the  quantity  o.^/Ab  where 

Oo  is  a  representative  size  of  a  turbulent  element,  X  is  the  wave 

length  of  light,  and  r  is  the  distance  from  the  turbulence  to  the  observer. 

If  »  I  the  geometrical  optics  approximation  is  valid. ^ 

For  a  slit  of  width  "a"  the  size  of  the  diffraction  pattern  is  equal  to 
the  geometric  image  of  the  slit  when  q^/  A  ^  =1.  Thus  the  criterion 

above  is  seen  to  be  plausible. 

The  problem  of  transmission  through  an  inhomogeneous  medium  has  been 

treated  by  (using  the  geometrical  optics  approximation)  Bergmann^  and 
5 

Chandrasekhar  .  The  results  of  this  paper  are  consistent  with  some  of  those 
of  Chandrasekhar  but  are  derived  from  less  restrictive  assumptions. 

For  an  investigator  interested  in  image  transmission  the  distortion 
property  of  an  inhomogeneous  medium  is  usefully  described  in  terms  of  the 
distortion  of  a  point  source.  More  precisely,  let  a  given  point  in  the  object 
space  be  designated  by  the  pair  of  direction  components 

3  Booker,  H.  G.,  and  Gordon,  W.  E. ,  "Outline  of  a  theory  of  radio  scattering 
in  the  troposphere,"  J.  Geophys.  Res.  241  (1950). 

^  Bergmann,  P.  G.,  "Propogation  of  radiation  in  a  medium  with  random  inhomo¬ 
geneities,  "Phys.  Rev.  70,  486  (1946), 

5  Chandrasekhar,  S.,  "A  statistical  basis  for  the  theory  of  stellar  scintilla¬ 
tion,"  Monthly  Not.  Roy.  Astron.  Soc.  U2,  475  (1952). 


Reduction  of  Contrast 


-  3  - 


and  let  the  point  be  assigned  a  unit  isotropic  radiance*  The  exact  form  of 

4*3)  will  be  given  later.  If  the  distortion  action  of  the  medium 

changes  sufficiently  rapidly  over  an  interval  of  time,  the  successive  images 
fl'*,')',)  of  appear  essentially  as  a  continuum  of  points 

(a  blur)  in  the  image  space.  For  the  purposes  of  this  discussion,  a  function 
CPis  defined  which  assigns  to  each  object  point  (  4^3  )  the  time-averaged 

radiance j  5^3^  of  the  image  point  .  The  time-averaged 

radiance  distribution  N  on  the  image  space  is  then  representable  as  a  con¬ 
volution  of  the  function 6^  with  the  radiance  distribution  N  on  the  object 
space: 


^ 


(1) 


where  is  the  object  space,  namely  the  collection  of  all  directions 

.  Since  the  location  of  the  observer  is  held  fixed  throughout  the 
entire  discussion,  the  location  symbol  x  has  been  suppressed  in  the  radiance 
function  notation  (ordinarily  written  Ni^,  ) 

In  the  following  discussion  attention  is  given  tc  the  statistics  of  the 
deviations  undergone  by  a  ray  traversing  a  path  betvreen  a  point  source  and  a 
screen  some  distance  away.  Because  of  an  optical  reciprocity  law  this 
derivation  applies  equally  well  to  an  extended  field  of  view,  from  which  light 
is  received  by  a  'point'  receiver  such  as  an  eye  or  a  lens. 

IffiTHOD  OF  APPROACH 

A  theoretical  treatment  for  transmission  through  any  inhomogeneous 


Reduction  of  Contrast 


-  4  - 


medium^s-  has  been  evolved  which  is  based  on  the  followj.ng  considerations i 
The  path  of  sight  is  assumed  to  lie  within  a  (spatially  and  temporally) 
uniformly  turbulent  medium  made  up  of  uncorrelated  optically  inhomogeneous 
regions  with  dimensions  small  compared  to  the  distance  through  which  the 
ray  travels.  Each  light  ray  undergoes  a  large  number  of  deflections.  Only 
small  angle  deflections  need  be  considered.  Using  only  these  assumptions  we 
may  consider  the  following-. 

Figure  1  shov;s  the  path  of  a  light  ray  emitted  at  traversing  a  turb¬ 
ulent  medium  toward  an  arbitrarily  oriented  screen  TT'  a  distance  r  from  the 
source.  The  plane  of  the  figure  is  perpendicular  to  TT'  and  passes  through 
the  path  P^Q  of  the  undeviated  ray.  Two  perpendicular  axes — the  x-  and 
y~axis  lie  in  TT'.  A  segment  TQT'  of  the  x-axis  is  shown.  The  y-axis  is 
normal  to  the  plane  of  the  figure  and  goes  through  Q,  The  ray  now  travels 
from  Pq  to  a  point  P^  determined  by  the  condition  that  the  ray  will  have 
turned  through  an  angle  whose  projection  on  the  plane  of  the  figure  is  of  a 
amall  fixed  magnitude  oo  (in  radians).  In  accordance  with  our  assumptions, 
the  distance  along  the  segment  PqP]_  is  large  compared  to  the  dimensions  of 
tlie  centers  of  turbulence.  The  projection  (in  the  manner  shown  in  the  figure) 
of  P-j_  falls  on  TT'  a  distance  x^  away  from  Q,  As  the  ray  continues  along  its 
path  from  P^^  it  ;«.ll  undergo  further  deflections,  '’/hen  its  direction  has 
been  changed  through  another  angle  whose  projection  in  the  plane  of  the 
figure  is  of  the  same  magnitude  o<  ,  P^  is  projected  on  the  screen,  where 
it  falls  a  distance  x^  away  from  the  projection  of  P^. 

•*!'Including  the  oceans  and  other  bodies  of  water,  wherein  refractive  inhorao- 
genities  may  result  from  suspensoids,  including  transparent  planktonic 
organisms,  in  addition  to  any  inhcmogenities  of  the  water  itself. 
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In  the  same  way  the  remainder  of  the  path  of  the  ray  is  divided  up  into 
segments  representing  equal  angles  of  deviation  oo  ,  and  the  projection  of 
each  endpoint  is  mapped  onto  the  screen  TT'.  Finally  will  be  the  x- 
distance  on  TT'  of  the  projection  of  from  the  projection  of  P^  where 

is  the  last  such  path  segment  traversed  by  the  ray  before  it  strikes 
the  screen.  No  generalit3'-  is  lost  if  P^  is  taken  to  lie  on  TT' ;  Since  the 
distance  over  large  compared  to  the  dimensions  of  the  centers 

of  turbulence,  it  follows  that  the  curvature  of  the  ray  between  P^^  and  P^ 
is  uncorrelated  with  its  curvature  during  its  past  or  future  history  so  that 
the  directions  of  projected  path  segments  are  uncorrelated. 

The  magnitude  of  x^  can  be  found  in  the  following  way,  using  the  small 
angle  assumption.  If  is  tlie  distance  from  P^  to  TT'  ,  then  x^  =  ,  For 

a  large  number  n  of  deflections,  the  value  of  will  approach  an  average 

value  given  by  (n  -  i)  d^,  vfhere  d^  is  the  average  distance  between  P^  and 
^  ^  .  Then  Xjj^  =  (n-i)ocd^.  The  algebraic  sign  associated  with  any  parti¬ 
cular  Xj^  has  equal  probability  of  being  +  or  It  can  then  be  seen  that 
the  distance  from  the  point  Q  to  the  point  Pj^  can  be  found  by  summation  over 
all  the  x^'s  vdth  the  appropriate  signs. 

Because  the  Xj^'s  are  uncorrelated,  they  can  be  considered  to  be  the  steps 
in  a  one  dimensional  random  walk:  Let  x^j,  the  last  of  the  x^'s  be  the  (signed) 
length  of  the  first  step  in  the  random  walk,  the  length  of  the  second  and 

so  on.  Therefore  absolute  magnitude  of  the  jth  step  will  be  -i  V=/  fh  • 

The  size  of  the  steps  in  this  one  dimensional  random  walk  thus  increase 
linearly  with  each  step. 
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TtiE  DETAILS  OF  THE  DERIVATION 


The  probability  distribution  for  a  one  dimensional  random  walk  ca.n 
be  calculated  in  the  following  way..  The  mean  square  deviation  Cr^(ol 
of  the  distribution  for  a.  walk  of  n  steps  is  equal  to  the  sum  of  the  mean 
square  deviations  for  each  step: 

<r‘(r<)  =  Zj.,  <rj  > 

where 


is  the  (signed)  length  of  the  jth  step. 
Therefore 

CT'Cnd  —  X.]L, 


n(n+i)  (an-v-O 

€ 


For  large  n  the  expression  in  brackets  approaches  :  0^/3  * 

Therefore 

(T^Cn)  =  ^  Y3  ^ 

6 

From  a  special  central  liiiiit  theorem  due  to  Liapounoff ,  it  can  be 
shovm  that  the  angular  probability  dj. stribution  approaches  a  normal  distri¬ 
bution  whose  mean  is  the  position  on  the  screen  reached  by  the  undeviated 
VB.J.  Liapounoff's  Theorem  states  the  follovdng:  let  x^,  Xp,  x^,...,  be 
independent  i-andom  variables,  and  denote  nfij  and  Oj  as  the  mean  and  standard 
deviation  of  x.^.  Suppose  that  the  third  absolute  moment  of  x.;  about  its 
raean: 

=  Expectation  (  I  m  j  I 

6  Cramer,  H.,  Mathematical  liethods  of  Statistics  (Princeton  Univ.  Press, 
I9Z1.6),  p  p.  215-216 , 
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is  finite  for  every  j  and  define  as; 

Then  if  the  condition  satisfied,  the  result  is  that 

the  sum  -oc  =  Xj  is  asymptotically  noirnal. 

This  can  be  applied  to  the  case  of  the  above  random  wallc  in  the  follow¬ 
ing  way: 

cn)  =  0^3  J  =  d.,,  [ - ^ -  J  . 


For  large  n: 


I A  ^ 


and, 

and  thus  the  hypothesis  of  Liapounoff's  theorem  is  satisfied. 

Hence  it  is  sho^^m  that  a  contiuum  of  rays  starting  vrlth  a  given  initial 
direction  and  which  traverse  a  medium  whose  optical  properties  are  uniformly 
turbulent  will  give  rise  to  a  continuum  of  terminal  points  on  the  target 
plane  such  that  the  density  of  the  points  projected  onto  the  x-axis  follows 
a  normal  distribution. 

Since 

I-  =  nd,,,  J 

vre  may  i/rite  alternatively 


for  fixed  . 
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V  3  ol. 


Thus  the  spread  of  the  Gaussian  probability  dist  ibution  increases  as  the 
3/2  power  of  r.  VJe  now  introduce  the  angular  variable  x/r.  In  terms 

of  ^  ,  the  spread  of  the  probability  distribution  Pj,  (S^)  is  proportional 
to  h .  The  normalized  Gaussian  probability  distribution  with  root  mean 
square  deviation  given  by 

CTxCh) 
is: 

The  above  derivation  may  be  extended  to  the  two-dimensional  case.  The 
new  derivation  is  similar  in  all  details  to  the  former.  In  particular  the 
same  angle  oC  is  used.  However,  the  projections  of  all  ray  segments  are  now 
onto  the  y-axis  in  the  plane  TT',  and  in  general  a  new  average  distance  dy 
must  be  introduced.  The  result  is: 


P.c'k)  - 


In  general  the  two-dimensional  probability  distribution  is  of  the  form: 

>)>«  po'fcf. 


r  -I  - r f ii_y_ „ r 

OiCHO-yG-l  lOi 


GiC'P=‘,  “Pa)  3ir 

where  ^  is  the  correlation  coefficient  between  Pj.  (  */x)  ^ 

Along  a  horizontal  path  of  sight  the  vertical  displacement  of  a  point  would 


be  expected  to  be  uncorrelated  vdth  the  horizontal  displacement.  In  this 
case 

I  (  I 

(Ph  (  )  =  PTT  (J^a)  CTyCH  1  ^ 

If  the  turbulence  is  isotropic,  i.e.,  if  d^^  =  dy  =  d,  then 


va^Cf))  J 


O-xCO  ^(T^Ch)  = 


_ h 

3  d 


o-ct)  . 
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If  we  set  ,  then  under  these  assumptions  (^=  o, 

we  have 

I  (  I  f  '!' 


air  <r^cH 


J. 


The  foregoing  fomulas  give  the  probability  that  the  terminal  point  on 


TT'  of  a  light  ray  aimed  along  P^Q  lies  in  the  direction  ^  seen 

from  Pq.  By  tracing  the  light  ray  in  the  opposite  direction,  these  formulas 
admit  the  folloi'dng  interpretation:  they  give  the  probability  that  the 
image  of  the  point  on  TT'  is  at  the  origin  (o,o)  of  TT'*  For  fixed 

r,<?r  is  thus  the  required  function  defined  in  (l), 

OPTICAL  AIR  STATE 

The  lovrer  limit  on  the  size  of  o<:  is  determined  by  the  condition  that 
the  direction  of  successive  steps  be  uncorrelated.  This  condition  will  be 
fulfilled  if  (as  assumed  above)  the  average  distances  d^^.  and  d  between 

y 

successive  Pj_'s  is  large  compared  to  the  dimensions  of  the  inhomogeneities 
of  the  medium.  It  follows  that  the  quantities  A^  =  A  jr  ~  ('<V3dy.) 

are  properties  only  of  the  medium  and  not  of  the  target.  The  pair  of  num¬ 
bers  (A„,A„)  is  called  the  optical  air  state  of  the  medium  and  gives  a  useful 
y 

measure  of  the  optical  turbulence  of  the  atmosphere.  If  the  turbulence  is 
isotropic,  the  optical  air  state  is  a  single  number  A  =  A'x  =  Ay.  In  most 
applications  the  turbulence  is  essentially  isotropic.  The  assumptions^  =  o. 
Ay.  =  Ay  v/ill  be  in  force  for  the  remainder  of  the  paper.  As  a  consequence 
we  may  \\'rite 


p,(A)  = 


p.  = 


\/2Tr  A  r" 


ex. 


]  ■> 


exp 


I  J 
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and 

where 

—  C£>  <  ‘fi)  <  +  CX?  j 

TIlfE-AVERAGED  APPARENT  RADIANCE 

Actual  computations  of  the  time-averaged  apparent  radiance  distribution 
Nj,  can  now  be  carried  out  by  means  of  formula  (l)  in  which  the  explicit 
expression  (2)  for  the  function(p  is  used.  A  particularly  useful  formula 
for  Np  arises  in  the  following  special  but  important  context:  let  a  target 
occupy  a  subregion  of  object  space  '■  —  <=t<<  ooj  -oo<^ 

let  the  apparent  radiance  distribution  N^  over  ]Pjl  be  uniform  and  of  magni¬ 
tude  ^Nj..  Let  the  remaining  portion  of  the  object  space  (the  background 

of  the  target)  have  uniform  apparent  radiance  of  magnitude  Then, 

from  (1) : 

=  i r  ^1- =  I.  +  I  b^/^-  0^  -  L  t A'h (?i- 

-f  %  ■’fy  ti  u 

This  result,  in  more  detail,  states  that 

for  all  in 

A  perusal  of  the  derivation  of(?r  would  show  that  (3)  (and  in  fact  any 
formula  making  use  of  the  derived(S?j,)  is  strictly  applicable  only  to  those 
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regions  of  which  are  angularly  small.  It  is  precisely  this  kind 

of  target,  however,  which  is  of  primary  concern  in  contrast  reduction  problems 
involving  atmospheric  boil.  These  targets  have  their  time-averaged  radiance 
greatly  influenced  by  their  background  radiance.  On  the  other  hand,  the 
time-averaged  background  radiance  of  such  targets  is  virtually  unaffected 
by  the  target's  radiance. 

TIlffi-AVERAGED  APPAREffl’  CONTRAST 


The  introduction  of  the  notion  of  contrast  greatly  facilitates  the 
solution  of  visibility  problems.  Let  ]  and  be  any 
two  points  in  either  image  or  object  space.  The  contrast  q  (  ) 
of  (  )  with  respect  to  (  ‘/'i,  is  defined  as 


v/here  N  is  a  given  radiance  distribution,  >-  O » 

As  an  example  consider  the  apparent  radiance  distribution  defined  in  the 
preceding  section.  Let  (^x,  be  an  arbitrary  fixed  point  in  Then 

the  apparent  contrast  function  Cj.  on  has  the  following  form: 


Cl. 

O 


foh  ia  fi  , 

fob  (fc,  «  %  . 


In  this  context  we  can  say  that  the  target  has  a  contrast  Cr  with 

respect  to  its  background  In  accordance  with  the  general  definition 

of  contrast,  the  time-averaged  apparent  contrast  function  Cr  necessarily 
has  the  form: 
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Now  let  be  an  arbitrary  fixed  point  in  .  In  view  of  the 

comments  of  the  preceding  section,  lire  may  for  all  important  cases  write 
CyL'^x,%)  as: 

where  is  any  point  in  ^  .By  adopting  this  convention  (3)  'i~ 

yields  the  folloxiring  approximate  but  highly  useful  contrast  formula: 

%  ^5) 

From  this  we  immediately  deduce  the  g  neral  facts  that 

^  1C^\  forall  int  , 

and  the  handy  rule  of  thumb  for  very  small  targets  of  area  Q.  and  arbitrary 
shape: 

^  ^  a 

Cp-  •  (6) 

This  points  up  among  other  things  an  inverse-cube  decay  of  Cj.  with  range  r. 
This  relation  may  be  used  to  explain  in  simple  terras  the  loss  of  fine  detail 
on  targets  as  the  observer- target  distance  is  increased.  For  by  subdividing 
the  given  target  into  sufficiently  small  elements  each  of  area  ^2-  ,  (6) 

ivill  describe  the  time-averaged  contrast  of  the  element  against  its  back¬ 
ground  consisting  of  the  remainder  of  the  target  area.  Now  while  the  flux 
from  the  element  falls  off  at  a  rate  proportional  to  ,  the  contrast 

(and  hence  the  detectability)  of  the  element  falls  off  at  a  rate  proportional 
to  .  In  short,  fine  detail  on  a  receding  target  is  lost  at  a  greater 

rate  than  the  flux  from  the  entire  target.  As  a  consequence,  the  detail  on 
a  target  is  lost  from  view  sooner  than  the  target  itself. 
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Despite  the  concession  to  comp?.exity  made  in  (4)j  the  resulting  formula  (5) 
faithfully  reproduces  the  main  features  of  the  time-averaging  affect  of 
atmospheric  boil  phenomena,  as  the  following  examples  illustrate. 


EXAilPLES 


A.  Long  Straight  Edge.  Consider  an  infinitely  long  straight  boundary  between 
two  regions  of  object  space  each  mth  a  distinct  uniform  apparent  radiance 
(Figure  2).  For  this  example,  the  target  consists  of  all  points 
( L,  'i',)  such  that  ^  0  i  The  background  is  the  complementary 

region  of  all  points  (*/*■,  )  such  that  ‘/(x  >  O  .  is  assigned  the 

uniform  apparent  radiance  ^ .  From  (3)  ivith  (2); 

<5  a? 


CO 

=  Lw,-  /  exf/-  ^  j  dPL  J  *  tVf  • 


M_  is  plotted  in  Figure  2. 


The  time  averaged  apparent  contrast  is,  by  (5): 

The  limits 


GJCp 

V<^..  r^r  ^ 


/3 

>c 


Ck  =  Cf.  ^ 

A^/hich  hold  for  all  ,  agree  with  the  values  one  would  expect  on  the  basis 

of  the  values  of  the  unshimmered  contrast  Q-r  discussed  in  the  preceding  action. 
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B.  Long  Thin  Bar.  A  long  thin  bar  may  be  represented  by  the  sub- 
region  of  consisting  of  all  points  (i'x,  )  such  that  —  o.<  . 

All  other  points  of  iF"  form  the.  background  .  '$).  is  assigned  the 

uniform  radiance  i  N).  ;  ,  the  radiance  b h  •  Then  from  (5)  with  (2): 


a.  o« 


It  is  easy  to  verify  that  j  C^j  attains  a  maximum  at  the  point  (0,0), 

,C^(0,0)  is  the  extreme  contrast  attained  by  Cj.  over  Y  >  i.e,,  Cj.  (0,0)  is 
either  a  maximum  (if  Cr(0,0)  is  positive)  or  a  minimum  (if  Cj.(0,0)  is  negative). 
It  is  clear  from  the  geometry  in  this  and  the  preceding  example  that  C^  is 
independent  of  ,  Some  contrast  profiles  of  shimmered  bars  are  shown  in 
Figure  5.  Observe  that  the  widths  of  the  bars  are  given  in  terms  of  the 
dimensionless  parameter  .  These  facts  will  be  of  some  use  in  the 

discussion  of  the  measurement  of  A, 

C.  Rectangle.  A  rectangle  in  centered  on  the  origin  is  defined 

as  the  subregion  consisting  of  all  points  such  that  <x  , 

•  From  (5)  mth  (2),  after  some  reductions. 


The  computations  of  Cj.  for  a  rectangle  can  therefore  be  carried  out  by 
using  the  results  of  exajnple  B, 
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D.  Circular  Disk.  A  circular  disk  in  centered  on  the  origin 

and  of  angular  radius  ^  is  the  set  of  all  points  such  that 

(Figure  3).  Prom  (5)  with  (2)  the  associated  contrast  expression  is: 

An  examination  of  the  function  Gj.  would  show  that  in  this  case  C^.  j  o}— 

Cyi^sc,  (circular  s3mimetry) .  It  is  therefore  sufficient  to 

consider  the  evaluation  of 


Ct-  <  5^=rj  —  STCA  a 


(Px-K)^-4- 

A 


(2 


Adopting  the  follov;ing  transformation  of  coordinates:  —  ’P'cosO'^  ^ 


we  have: 


Ci- 


3TCf{  h 


exp 


The  inner  integral  is  readily  evaluated: 
Sir  IT 


c'oso-' 


S/0&')  ~ 


/^h  ? 


where  lo  is  a  modified  Bessel  function  of  the  first  kind,  of  zero  order,  which 
has  the  series  representation: 


loC^)^ 


oo 


ill 

(A!)^ 
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Hence 


C^('/(.,o)==  ^  J  e^cpj-  2/i  /- 


which  may  be  rewritten  as: 


c.iKo)  =  ijrp 


CL 


whehe  PC^-;  ^+i)=  j"  defines  the  incomplete  Gamma  function,  for 

which  tabulations  exist^,  and  4'^—  4  -  4^/2.  A  h  .  It  is  easy 

to  check  that  1  Cj.  j  is  a  decreasing  function  for  increasing  s4c  >  so  that 
j I  attains  a  maximum  at  (0,0) i  Thus  if  jc^  (0,0) |  is  below  a  threshold 
contrast,  the  entire  disk  is  undetectable.  The  formula  for  Cj.(0,0)  is 
readily  obtained: 

_  h 

C^(0;0)=:  r(^y2Ah  i  l)  =  C^(l~e  ). 


a 


For  angularly  small  disks 

C^CO.o)  —  Ch  2-rr/\h3  ’ 

where  GL=~rr(4h)^  ,  and  which  agrees  with  (6).  If  ,  i.e.,  if  the 

disk  is  infinitely  large,  then  since  p(co  jk-^\)=:Jk  !  >  "'^he  general  expression 

yields: 

-  9^  __  ou 


C,e 


4  =  0 


-  =  Cy  ^ 

(A!)  ’ 


as  expected. 

E'.  One  example  of  the  effect  of  atmospheric  boil  on  the  ability  to  resolve 
two  objects  which  are  close  together  is  illustrated  in  Figure  4.  Two  identical 
black  bars  of  angular  -vidth  /A  h  ,  sepai-ated  an  angular  distance  , 

have  their  contrast  profiles  smeared  into  a  characteristic  shape 


’  Pearson,  Karl,  Tables  of  the  Incomplete  Gamma  function  (Camb.  Univ.  Press, 

1946) . 
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shown  by  the  dashed  curve.  The  object  space  contrast  of  -1  at  the  center  of 
each  bar  is  invariably  raised  to  -  0.6  in  the  image  space,  and  the  zero 
contrast  at  the  center  of  the  gap  between  the  bars  is  invariably  pulled  down 
to  -  0.74.  Such  predictable  features  of  shimmered  contrast  profiles  can  be 
used  to  generate  an  experimental  procedure  for  detemining  A,  one  of  which  is 
appended  belowi 

MEASURMNT  OF  OPTICAL  AIR  STATE 


For  a  bar  whose  width  2a  is  less  than  5\[/\^  >  the  values  | 

are  measurably  less  than  |  .Cr .(.  Figure  5  depicts  the  contrast  profiles 
of  four  long  black  bars  (hence  Cr  =  -1)-.  The  width  of  each  differs 
by  a  factor  of  four  from  its  predecessor.  From  example  B, 

cu 

-  I 


By  means  of  this  integral,  each  choice  of  the  ratio  a<x/\f/C\  j-  leads  to  a 
corresponding  value  0^(0, 0).  These  correspondences  are  graphed  in  Figure  6. 

The  experiraental  procedure  consists  in  measuring  Cj;(0,0)  for  each  bar  by 
scanning  the  bar  with  a  telephotoraeter  at  a  given  range  r.  The  four  corres¬ 
ponding  values  of  are  then  picked  from  the  graph.  Since  r  is 

fixed,  and  2a  is  Icno^vn  for  each  bar,  A  can  then  be  determined, 

’'le  observe  finally  that  before  using  the  graph  in  Figure  6  the  measured 
values  Cj.(0,0)  must  be  corrected  to  remove  the  reduction  effect  on  the  contrast 
induced  by  the  scattering  and  absorption  mechanisms  of  the  atmosphere.  For 
horizontal  paths  of  sight,  with  the  usual  uniformity  properties,  the  corrected, 
values  are  given  by  C^0,0)  exp^-h  t/ij  v/here  L  is  the  attenuation  length 
of  the  atmosphere. 
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CAPTIONS 


Illustrating  the  geometry  of  the  derivation  of  .-(P. 

Illustrating  the  time-averaged  radiance  profile  of  a  shimmered 
edge  (Example  A). 

Diagram  for  the  derivation  in  example  D , 

The  characteristic  shape  of  the  time-averagedcontrast  profile 
of  two  long  black  bars  of  the  dimensions  shown. 

Illustrating  the  characteristic  time-averaged  contrast  profiles 
of  long  black  bars  as  a  function  of  bar-width. 

Graph  used  in  the  experimental  determination  of  A. 


Figure  6 
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